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ABSTRACT 

Galaxy source counts in the infrared provide strong constraints on the evolution of the bolometric energy 
output from distant galaxy populations. We present the results from deep 24//m imaging from Spitzer surveys, 
which include « 5 x lO"* sources to an 80% completeness of ~ 60 /iJy. The 24^m counts rapidly rise at near- 
Euclidean rates down to 5 mJy, increase with a super-Euclidean rate between 0.4-4 mJy, and converge below 
■^0.3 mJy. The 24 /im counts exceed expectations from non-evolving models by a factor > 10 at 5^ ^ 0. 1 mJy. 
The peak in the differential number counts corresponds to a population of faint sources that is not expected from 
predictions based on 15/im counts from ISO. We argue that this implies the existence of a previously undetected 
population of infrared-luminous galaxies at z ~ 1-3. Integrating the counts to 60 ^Jy, we derive a lower limit 
on the 24/im background intensity of 1 .9 ± 0.6 nW m~^ sr~' of which the majority {^ 60%) stems from sources 
fainter than 0.4 mJy. Extrapolating to fainter flux densities, sources below 60 /zJy contribute 0.8^q4 nW m~^ 
sr~' to the background, which provides an estimate of the total 24 /im background ofl.l^Q^ nW m~^ sr~' . 



Subject headings: cosmology: observations — galaxies: evolution 
photometry — infrared: galaxies 



galaxies: high-redshift — galaxies: 



1. INTRODUCTION 

From the first detections of infrared (IR) luminous galaxies, 
it was clear that they represent phenomena not prominent in 
optically selected galaxy surveys (e.g., Rieke & Low 1972; 
Soifer, Neugebauer, & Houck 1987). Locally, galaxies ra- 
diate most of their emission at UV and optical wavelengths, 
and only about one-third at IR wavelengths (5- 1000 /im; 
Soifer & Neugebauer 1991). IR number counts from ISO 
indicate that the IR-luminous sources have evolved rapidly, 
significantly faster than has been deduced from optical sur- 
veys, which implies that IR-luminous galaxies make a sub- 
stantial contribution to the cosmic star-formation rate density 
(e.g., Elbaz et al. 1999; Franceschini et al. 2001). The detec- 
tion of the cosmic background by COBE at IR wavelengths 
shows that the total far-IR emission of galaxies in the early 
Universe is greater than that at optical and UV wavelengths 
JFixsen et al. 1998; Hauseretal. 1998; Madau & Pozzett] 
llpOO; Franceschini et al. 2001|), which suggests that a large 
fraction of stars h ave form ed in IR-luminous phases of 
galaxy activity (Elb az et all 12002). Studies from ISO at 
15 /^m and 170 /im have inferred that the bulk of this back- 
ground originates in disc rete sources with z < 1 .2 ( Dole et al. 
1200 It lElbaz et all l2002h . However, the peak in the cos- 

' This work is based on observations made with the Spitzer Space Tele- 
scope, which is operated by the Jet Propulsion Laboratoiy, California Institute 
of Technology under NASA contract 1407. 

- Steward Observatory. The University of Arizona, 933 North Cherry Av- 
enue, Tucson, AZ 85721; Electronic Address: papovich@as.arizona.edu 

^ Institut d'Astrophysique Spatiale, bat 121, Universite Paris Sud, F91405 
Orsay Cedex, France 

* Departamento de Astrofisica Molecular e Infrarroja, IBM, CSIC. Serrano 
113b, 28006 Madrid, Spain 

"^ Michelson Science Center, California Institute of Technology, Pasadena, 

CA 91109 

^ Space Science Institute, 4750 Walnut, Suite 205, Boulder, CO 80301 
' National Optical Astronomy Observatory, 950 North Cherry Avenue, 

Tucson, AZ 85719 

* Department of Physics, University of California, Berkeley, CA 94720 



mic IR background extends to r^ 200/im dFixsen et alJll99« 
Hauser & Dwek 2001). The spectral-energy distributions 
(SEDs) of luminous IR galaxies (L[r ~ lO' '" L©) peak be- 
tween 50-80 /im ( Dale et al. 2001 ). If these objects consti- 
tute a major component of the cosmic IR background, then 
it follows there is a significant population of IR-luminous 
galaxies atz'--^ 1.5-3 — distances largely unexplored by ISO. 
The mid-IR 24 /im band on the Multiband Imaging Pho- 
tometer for Spitzer (MIPS, Rieke et al. 2004) is particu- 
larly well suited for studying distant IR-luminous galax- 
ies. Locally, the mid-IR emission from galaxies relates 
almost linearly with the total IR luminosity over a range 
of galaxy types (e.g., Soinoglio et al. 1995: Cha rv & Elbazl 
2001; Roussel et al.'2001; Papovich & Bell 2002), and there 
are in dications this holds at higher redshifts (e.g., Elbaz et ^ 
l2002h . Because the angular resolution of Spitzer is signifi- 
cantly higher for the 24 /im band relative to 70 and 160 /im, 
the 24 /im confusion limit lies at fainter flux densities. This 
allows us to probe the IR emission from many more sources 
and at higher redshifts than with the MI PS long er wavelength 
bands (e.g.. IPapovich & Beiill2002; .Dole et alJl2003ft . Here, 
we present the number counts of ~ 5 x lO'* sources detected 
at 24 /im in deep Spitzer surveys, and we suggest that the faint 
Spitzer detections probe a previously undetected population 
of very luminous galaxies at high redshifts. Where applica- 
ble, we assume 17„, — 0.3, f^A = 0.7, and Hq = 70 km s~' 
Mpc"^ 

2. THE data and SOURCE SAMPLES 

The data used in this work were obtained in five fields from 
early Spitzer characterization observations and from time al- 
located to the MIPS Guaranteed Time Observers (GTOs). Ta- 
ble ^ lists their properties and 24 /im source densities. The 
MIPS fields used here subtend the largest areas (~ 10.5 deg^) 
and widest range in flux density (50-0.06 mJy) available to 
date from the Spitzer mission, and therefore are the premier 
dataset for studying the mid-IR source counts. 
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TABLE 1 
Properties of Deep, Spitzer Fields 





R.A. 


Decl. 


{!.) 


Area 


(fexp) 


Ctio% 


A'(>C8o%) 


N{> 300 ^Jy) 


Field 


(J2000.0) 


(J2000.0) 


(MJysr') 


(arcmin^) 


(sec) 


(MJy) 


(arcmin^- ) 


(arcmin ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Marano 


03 13 52 


-55 15 23 


19.7 


1296 


236 


170 


2.0 


0.9 


CDF-S 


03 32 28 


-27 48 30 


22.5 


2092 


1378 


83 


4.5 


0.7 


EGS 


14 16 00 


+52 48 50 


19.5 


1466 


450 


110 


3.4 


0.7 


Bootes 


14 32 06 


+34 16 48 


22.7 


32457 


87 


270 


1.0 


0.8 


ELAIS 


16 09 52 


+54 55 00 


18.2 


130 


3232 


61 


5.7 


0.6 



Note. — Col. (1) Field name; (2) right ascension in units of hours, minutes, and seconds; (3) declination in units of degrees, arcminutes, and arcseconds; (4) 
mean 24 fj,m background; (5) areal coverage; (6) mean exposure time; (7) 80% completeness limit; (8) source density with Si^ > C8o% . (9) source density with 
5^ > 300 ^Jy. 



The MIPS 24 /im ima ges were p rocessed with the MIPS 
GTO data analysis tool ^Gordon et aL.2004) . The measured 
count rates are corrected for dark current, cosmic-rays, and 
flux nonlinearities, and then divided by flat-fields for each 
MIPS scan-mirror position. Images are then corrected for 
geometric distortion and mosaicked. The final mosaics have 
a pixel scale of ~ 1"25 pix~', with a point-spread function 
(PSF) full-width at half maximum (FWHM) of ~ 6". 

We performed source detection and photometry using a set 
of tools and simulations. Briefly, we first subtract the im- 
age background using a median filter roughly four times the 
size of source apertures (see below), and use a version of the 
DAOPHOT software (Stetson 1987) to detect point sources. 
We filter each image with a Gaussian with a FWHM that is 
equal to that of the MIPS 24 /im PSF, and identify positive 
features in 10"-diameter apertures above some noise thresh- 
old. We then construct empirical PSFs using 20-30 bright 
sources in each image, and optimally measure photometry by 
simultaneously fitting the empirical PSFs to all sources within 
~ 20" of nearby object centroids. The source photometry cor- 
responds to the flux of these PSF-apertures within a diameter 
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Fig. 1. — Reliability of source detection and photometry for the CDF-S, 
one of the deep Spitzer fields. The thick line shows the fraction of artificial 
sources detected in the image as a function of input flux density. For the 
CDF-S this fraction is 80% complete at Cgo'H = 83 /xJy (indicated by the 
dotted line). The red line shows the ratio of the number of sources detected 
in a 'negative' of the 24fiTn image to the number of 'positive' sources as a 
function of flux density, which provides an estimate of the number of spurious 
sources due to noise features. 



of 37''4, and we apply a multiplicative correction of 1.14 to 
account for light lost outside these apertures. 

To estimate completeness and photometric reliability in 
the 24 /im source catalogs, we repeatedly inserted artifi- 
cial sources into each image with a flux-distribution approx- 
imately matching the measured number counts. We then re- 
peated the source detection and photometry process, and com- 
pared the resulting photometry to the input values. In fig- 
ure [2 we show for the Chandra Deep Field South (CDF-S; 
one of the deep, Spitzer fields) the relative fraction of sources 
within r <2'.'3 (half the FWHM) and with a flux difference 
< 50% compared to their input values. From these simu- 
lations, we estimated the flux-density limit where 80% of 
the input sources are recovered with this photometric accu- 
racy, and these are listed in table ^ Simultaneously, we es- 
timate that down to the 80% completeness limit the number 
of sources that result from fainter sources either by photomet- 
ric errors or the merging of real sources is < 10%. We also 
repeated the source detection and photometry process on the 
negative of each MIPS 24 /im image. This test provides an 
estimate for the number of spurious sources arising from the 
noise properties of the image, as shown in figure ^ For all 
of our fields, the spurious-source fraction for flux densities 
greater than the 80% completeness limit is < 10%. 

3. THE 24 ^m SOURCE COUNTS 

Figure|2shows the 24 /xm cumulative and differential num- 
ber counts that have been averaged over the fields listed in 
table n The differential and cumulative counts (corrected 
and uncorrected) are listed in table |2 The faintest datum 
(denoted by the open symbol in figure |2} is derived to the 
50% completeness limit for the European Large-Area ISO 
Survey (ELAIS) field (C5o% = 35 /iJy). The remaining (less 
deep) fields are used after correcting for completeness to the 
80% level only. Error bars in the figure correspond to Pois- 
sonian uncertainties and an estimate for cosmic variance us- 
ing the standard deviation of counts between the fields. For 
Si, < 83 /iJy where the counts are derived solely from the 
smaller ELAIS field, we estimate the uncertainty (18%) us- 
ing the standard deviation of counts at faint flux densities in 
cells of 130 arcmin^ from the CDF-S (see below). We have 
ignored the contribution to the number counts from stars at 
24 /im, which are negligible at these Galactic latitudes and 
flux densities based on preliminary Spitzer observations. 

At bright flux densities. Si, > 5 mJy, the differential 
24 /im source counts increase at approximately the Euclidean 
rate, dN /dSu ^ S~^^^, which extends the trends observed by 
the IRAS 25 /im population by two orders of magnitude 
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Fig. 2. — Cumulative (Left panel) and differential (Right panel) 24/im number counts. The differential counts have been normalized to a Euclidean slope, 
dN/dSi, ~ 5^"'^. The solid stars show the average counts from all the Spitz.er fields (see TablefTl. and corrected for completeness to their respective 80% limits. 
The open star corresponds to counts brighter than the 50% completeness limit from the ELAIS field. The en'or bars coiTespond to counting uncertainties and 
a cosmic-variance estimate based on the standard deviation of the field-to-field counts from the different fields. Each flux bin is A(logS,y) = 0.15 dex. The 
shaded diamonds coirespond to IRAS 25 fj,m number counts from Hackins & Soifer 1 1991, and adjusted assuming U24S^{24fj.m) = U25Si,{25fini)). The curves 
show the predictions from various contemporary models from the literature (see figure inset; and adjusted slightly to match the observed IRAS counts), and a 
model based on the local ISO 15^m luminosity function and assuming non-evolving galaxy SEDs. 



dHacking & Soife r'lQQl: "Shu pe et alJfTOQSl) . For S^ ~ 0.4- 
4 mJy, the 24 /xm counts increase at super-Euclidean rates, 
and peak near 0.2-0.4 mJy. This observation is simila r to the 
tren d observed in the ISO 15 /xm source counts ( Elbaz et al.l 
1199 91 but the peak in the 24 /xm differential source counts 
occurs at fluxes fainter by a factor « 2.0. The peak lies above 
the 80% completeness limit for nearly all fields, and is seen 
in the counts of the fields individually. Thus the observed 
turn over is quite robust. The counts converge rapidly with a 
sub-Euclidean rate at < 0.2 mJy. We broadly fit the faint-end 
(Sij ~ 35 - 130 /iJy) of the differential number counts using a 
powerlaw, dN/dS^ ^ C {S ^ / 1 mJy)'"' , with a =-1.5 ±0.1. 
This result is consis tent with a separate analysis based solely 
on the ELAIS field ( Chary et al..2004.) . 

The observed counts are strongly inconsistent with expec- 
tations from non-evolving models of the local IR-luminous 
population. In figure |2l we show the 24 /im counts derived 
from the local luminosity function at ISO 15 /im (Xu 2000) 
and assuming local galaxy SEDs (.Dale et al.l2001.) . We have 
used the the local ISO 15 /zm luminosity function, because 
the A:-correction between rest-frame ISO 15 ^m and observed 
MIPS 24 ijm bands is minimized at higher redshifts (z ^ 0.6), 
which is more appropriate for the counts at fainter flux densi- 
ties. However, using the local IRAS 25 /im luminosity func- 
tion dShupe et al.lll99a) yields essentially identical results. 
While the non-evolving fiducial model is consistent with the 
observed 24 /im counts for 5^ > 20 mJy, it underpredicts the 
counts at S^y < 0.4 mJy by more than a factor of 10. 

Because the deep Spitzer fields are measured in many sight- 
lines with large solid angle, we can estimate the fluctuations 
in the number of sources as a function of sky area and flux 
density in smaller-sized fields. For example, in the CDF-S, 
which achieves an 80%-completeness limit of 83 /iJy over 
~ 0.6 deg^ (see table [2 and figure 0, we have computed 
the variance in the number of sources in solid angles of 100 
and 300 arcmin^. For the sources that span the peak in the 



differential counts (0.1-1 mJy), the fluctuation in the num- 
ber of sources in flux-density bins of 0.15 dex is roughly 
15% in these areas. This implies that small-sized fields suf- 
fer sizable field-to-field variation in the number of counts 
from the cosmic variance of source clustering. This effect 
is present in even larger fields: the number density of sources 
brighter than 0.3 mJy varies by ^ 10% in fields of ^ 0.5 deg^ 
(table 0, and is consistent with fluctuations expected from 
galaxy clustering on fields of this size at z '-^ 1, (scale lengths 
of 20-70 Mpc). The counts presented here average over fields 
from many sightlines and significantly larger areas. We con- 
servatively estimate that variations due to galaxy clustering 
correspond to uncertainties in the number counts of a less than 
a few percent in each flux bin. 

4. INTERPRETATION AND DISCUSSION 

The form of the observed 24 /im source counts differs 
strongly from predictions of various contemporary models 
(see figure |2ji. Four of the models are phenomenologi- 
cal in approach, which parameterize the evolution of IR- 
luminous galaxies in terms of density and luminosity to 
match observed counts from ISO, radio, sub -mm, and other 
datasets. S everal of these models (Charv & E lbad 1200 It 
iRing & Row an-Robinson 2003; Xu et al. 2003) show a rapid 
increase in the number of sources at super-Euclidean rates at 
relatively bright flux densities (Si, > 10 mJy) and peak near 
1 mJy. These models generally predict a redshift distribution 
for the MIPS 24 /im population that peaks near z ^ 1, based 
largely on expectations from the ISO populations, and they 
overpredict t he 24 /xm nu mber counts by factors of 2 - 3 at 
^ 1 mJy. ThelLaga che et al.. C2003.) model predicts a roughly 
Euclidean increase in the counts for S,y > 10 mJy. The shape 
of the counts in this model is similar to the observed distribu- 
tion, but it peaks at 5,^ ^ 1 mJy, at higher flux densities than 
the observed counts. This model predicts a redshift distribu- 
tion that peaks near z ^ 1, but tapers slowly with a significant 
population of IR-luminous galaxies out to z > 2 dDole et all 
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120031) . 

The model of Balland, Devriendt, & Silk (2003) is based 
on semi-analytical hierarchical models within the Press- 
Schecter formalism in which galaxies identified as 'interact- 
ing' are assigned IR-luminous galaxy SEDs. This model in- 
cludes additional physics in that the evolution of galaxies de- 
pends on their local environment and merger/interaction his- 
tories. Although this model predicts a near-Euclidean in- 
crease in the counts for 5^ > 10 mJy, the counts shift to 
sub-Euclidean rates at relatively bright flux densities. The 
semi-analytical formalism seems to not include important 
physics that are necessary to reproduce the excess of faint 
IR sources. This illustrates the need for large-area multi- 
wavelength studies of Spitzer sources to connect optical- and 
IR-selected sources at high redshift to understand the mecha- 
nisms that produce IR-luminous stages of galaxy evolution. 

The peak in the 24 ^um differential number counts occurs 
at fainter flux densities than predicted from the phenomeno- 
logical models based on the ISO results. This may suggest 
possibilities such as a steepening in the slope of the IR lu- 
minosity function with redshift, or evolution in the relation 
between the mid- and total IR. Phenomenological models 
which reproduce the IR background predict a faint-end slope 
of the IR luminosity function that should be quite shallow 
at high redshifts, with 'L*' luminosities that correspond to 
LiR > 10^ ' Lq for z > 1 (see Hauser & Dwek 2001 ). For most 
plausible IR luminosity functions, galaxies with L* lu minosi- 
ties dominate the integrated luminosity density. JSlbaz et al.l 
(|2002) observed that the redshift distribution of objects with 
these luminosities in deep ISO surveys spans z — 0.8-1.2, 
and that these objects constitute a large fraction of the total 
cosmic IR background. Therefore, it seems logical that ob- 
jects with these luminosities dominate 24 /im number counts 
at 0.1 -0.4 mJy, and it follows that their redshift distribution 
must lie at z ~ 1-3 (i.e., where this flux density corresponds 
to ~ 10''~^^ Lq, using empirical relations from Papovich & 
Bell 2002). Indeed, a similar conclusion is inferred based on 
a revised phenomenological model using the 24 /im number 
counts presented here (Laaache et al. 2004), and allowing for 
small changes in the mid-IR SEDs of IR-luminous galaxies. 
Examples of MIPS 24 /im sources at these redshifts and lumi- 



nosities hav e been readily identified in optical ancillary data 
iLe Floc'h e t al. 2004). We therefore attribute the the peak in 
the 24 /im differential number counts at fainter flux densities 
to a population of luminous IR galaxies at redshifts higher 
than explored by ISO. 

Integrating the differential source-count distribution pro- 
vides an estimate for their contribution to the cosmic IR back- 
ground at 24 /im, i.e, Ii, ~ j dN /dS^S^dS^. For sources 
brighter than 60 /iJy, we derive a lower limit on the total back- 
ground of vIy{2Aiiva) = 1.9 ±0.6 nW m~^ sr~^ Due to the 
steep nature of the source counts, most of this background 
emission results from galaxies with fainter apparent flux den- 
sities. We find that 2± 60% of the 24 /im background origi- 
nates in galaxies with Si, < 0.4 mJy, and therefore the galaxies 
responsible for the peak in the differential source counts also 
dominate the total background emission. Our result is consis- 
tent with the COBE DIRBE upper limit vl„{25p.m) < 7 nW 
m~^ sr~ ' inferred from fluctuations in the IR background 
jKashli nskv & Odenwald 2000; Hauser & Dwek 2001). As 
a further estimate on the total 24 /im background intensity, 
we have extrapolated the number counts for S^, < 60 /iJy 
using the fit to the faint-end slope of the 24 /im number 
counts in § 3. Under this assumption, we find that sources 
with Sy < 60 /iJy would contribute 0.8;!;q 4 nW m~^ sr~^ to 
the 24 /im background, which when summed with the above 
measurement yields an estimate of the total background of 

vlj)°'\24-fj,m) = 2.7!i;i nW m'^ sr"'. For this value, the 
sources detected in the deep Spitzer 24 /im surveys produce 
~ 70% of the total 24 /im background. 
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TABLE 2 
Measured Spitzer 24 ^m Number Counts 



Log5i, 


dN/dS^ 


S{dN/dS„) 


Log Si, 


N{>S^) 


S{N{> S^)) 


(mJy) 


(mJy"' sr') 


(mJy"' sr') 


(mJy) 


(sr-') 


(sr') 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


-1.475 


2.5(1.4) xlO' 


4.7(1.5) xlO** 


-1.550 


2.1 (1.1) xlO** 


1.2(0.60) xlO' 


-1.325 


1.5(1.2) xlO' 


2.7 (0.97) X 10** 


-1.400 


1.2(0.92) xlO** 


7.6(4.9) xlO* 


-1.175 


8.7 (6.9) X 10« 


1.4(1.4) xlO* 


-1.250 


9.6 (7.2) X 10' 


3.7 (0.68) X 10« 


-1.025 


5.5(5.2) xlO« 


5.8(1.6) xlO' 


-1.100 


6.1 (5.6) X 10' 


5.9(5.6) xlO' 


-0.875 


3.1(3.1) xlO« 


7.1(4.9) xlO« 


-0.950 


3.9 (3.9) X 10' 


4.0 (4.0) X 10' 


-0.725 


1.7(1.6) xlO* 


7.6(6.0) xlO^ 


-0.800 


2.9(2.5) xlO' 


5.5(5.0) xlO' 


-0.575 


8.3(7.5) xlO' 


1.5(0.66) xlO« 


-0.650 


1.8 (1.6) X 10' 


1.8(1.6) xlO^ 


-0.425 


3.2 X 10' 


5.7 X 10^ 


-0.500 


8.6x10'' 


1.2x10^ 


-0.275 


1.2 x 10' 


1.9x10^ 


-0.350 


4.4 X lO*" 


6.2x10" 


-0.125 


4.1 X lO*" 


9.6 X lO"* 


-0.200 


2.3 X 10*" 


3.8x10" 


0.025 


1.4 X lO*" 


3.9 X lO"* 


-0.050 


1.2x10* 


2.3 X 10" 


0.175 


5.8 X 10' 


2.1 X lO"* 


0.100 


7.1 xlO' 


1.7x10" 


0.325 


2.4 X 10^ 


1.1 xlO"* 


0.250 


4.1 xlO' 


1.4x10" 


0.475 


1.0 X 10^ 


6.5 X 10^ 


0.400 


2.3 X 10^ 


1.0x10" 


0.625 


3.4 X lO'* 


2.8x10^ 


0.550 


1.3x10^ 


7.1x10^ 


0.775 


1.5x10"* 


1.6 xlO^ 


0.700 


7.7x10'' 


5.5 X 10^ 


0.925 


6.0 X 10^ 


8.6 xlO^ 


0.850 


5.2x10'' 


7.4 xlO^ 


1.075 


2.7 X 10' 


4.9 X 10^ 


1.000 


3.4x10" 


6.0x10' 


1.225 


1.2 xlO^ 


2.7 X 10^ 


1.150 


2.0x10" 


4.0 xlO' 


1.375 


6.2 X 10^ 


1.7 xlO^ 


1.300 


1.3x10" 


2.8x10=' 


1.525 


2.2 X 10^ 


8.3x10' 


1.450 


7.3 X 10^ 


2.1 X 10' 


1.675 


1.1 xlO^ 


5.0x10' 


1.600 


4.4 xlO^ 


1.5x10' 



Note. — Col. (1) Flux density for differential number counts; (2) Corrected differential counts 
and (3) uncertainty; (4) Flux density for cumulative number counts; (5) Corrected cumulative 
counts and (6) uncertainty. Numbers in parentheses give the uncorrected values. 



